Electric control of ferromagnetism in Mn-doped semiconductor heterostructures 
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The interplay of tunneling transport and carrier-mediated ferromagnetism in narrow semicon- 
ductor multi-quantum well structures containing layers of GaMnAs is investigated within a self- 
consistent Green's function approach, accounting for disorder in the Mn-doped regions and un- 
wanted spin-flips at heterointerfaces on phenomenological ground. We find that the magnetization 
in GaMnAs layers can be controlled by an external electric bias. The underlying mechanism is 
identified as spin-selective hole tunneling in and out of the Mn-doped quantum wells, whereby the 
applied bias determines both hole population and spin polarization in these layers. In particular we 
predict that, near resonance, ferromagnetic order in the Mn doped quantum wells is destroyed. The 
interplay of both magnetic and transport properties combined with structural design potentially 
leads to several interrelated physical phenomena, such as dynamic spin filtering, electrical control 
of magnetization in individual magnetic layers, and, under specific bias conditions, to self-sustained 
current and magnetization oscillations (magnetic multi-stability). Relevance to recent experimental 
results is discussed. 
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I. INTRODUCTION 

Electric control of magnetism in nanostructures must 
be viewed as an important milestone on our road map 
for successful realization of spintronic devices. Although 
most of the operations in such devices ultimately should 
be based on spin-only processes, i.e., processes not asso- 
ciated with (highly dissipative) electric charge transport, 
to gain best benefits from such designs, spin must be 
manipulated both during the input, control, and read- 
out stage and eventually be coupled to charge. Several 
schemes achieving this goal have been explored both at 
the quantum and semi-classical level, such as the electric 
distortion of the orbital wave function of spin carriers in 
inhomogeneous (effective) magnetic fields electric g- 
tensor control 0, Q, or spin torque transfer. Here 
we explore, on theoretical grounds, the influence of an 
electric bias on the ferromagnetic state and feasibility of 
electric control of ferromagnetism in Gai_a;_j,Alj,Mnj;As 
multiple quantum wells. Structural design, including ef- 
fective potential profiling and doping to position emitter 
and collector quasi-Fermi levels, as well as tunneling is 
used to control hole density and spin polarization within 
the Mn doped layers. 

Dilute magnetic semiconductors (DMS) have been re- 
alized by doping of conventional ZnS-structured semi- 
conductors with elements providing open electronic d or 
/ shells. This has added yet another degree of freedom 
to the rich spectrum of physical phenomena in semi- 
conductors available for material design with potential 
for technological applications. A prominent example 
is bulk Gai_a;Mna;As where Mn on the Ga sites pro- 
vides both an open d-shell with a local magnetic mo- 
ment and a hole which may establish ferromagnetic or- 
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dering among the Mn d-electrons, a mechanism known 
as carrier-mediated ferromagnetism. 8-10] The prefer- 
entially anti-parallel alignment of the 3/2 spin of the 
mobile holes with the 5/2 spin of the localized Mn d- 
electrons promotes ferromagnetic ordering of the latter 
below a critical temperature of up to ~ 150 K. Theoreti- 
cal work has confirmed strong hybridization between the 
5d Mn and 3p electrons in the ground state. [ll| The ef- 
fective hole-concentration-dependent exchange field lifts 
the spin degeneracy of the holes' energy bands and thus 
goes hand in hand with hole spin polarization. The Mn 
ions sitting on Ga sites act as acceptors and are believed 
to give rise to acceptor levels which lie about 100 meV 
above the valence band edge.0,|^[l3| Photoluminescence 
experiments indicate the co-existence of holes bound to 
Mn sites and itinerant holes which participate in estab- 
lishing magnetic order amongst the Mn ions below T^.[l3| 

Since structural defects of bulk and confined layers of 
Gai_a;Mn2:As depend on grow conditions, Mn concentra- 
tion X , and annealing procedures it is not too surprising 
that experiments have come up with somewhat different 
conclusions regarding the "electronic structure of bulk 
Gai_a;Mn2;As" . More recent work seems to hint at the 
existence of an impurity band which forms at Mn concen- 
trations above 1.5 % leading to a metal-insulator tran- 
sition in high-quality GaMnAs. [T3 - [T6| The Fermi level 
in these samples is reported to lie in the impurity band 
and the valence-band properties remain largely GaAs- 
like.fl^ The radius of the Mn acceptor wave function has 
been measured to be about 2 nm, indicating that Mn^a 
is not a shallow acceptor, [l^ In contrast, other studies 
rather hint at a disordered top valence band edge con- 
taining the Fermi energy, but no isolated impurity band is 
present. Recent theoretical work has led to the conclu- 
sion that a tight -binding approach (within the coherent- 
potential approximation for disorder) and local-density 
functional theory -I- Hubbard U correction cannot ac- 
count for an isolated impurity band.fivj Other theoret- 
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ical work has lead to the conclusion that disorder may 
enhance ferromagnetic stability, [isl [T9| Ionized impurity 
scattering seems to play the dominant role in explaining 
Hall resistivity data.|20| 

Controlled growth of heterostructures containing crys- 
talline layers of GaMnAs of high structural quality has 
remained a challenge up to date. Nevertheless, tunnel- 
ing spectroscopy has confirmed siz e q uantization effects 
in GaMnAs quantum well lavers.|l6l [2ll [2^ However, 
compared to crystalline GaAs well layers, in an oth- 
erwise identical structure, the signature appears to be 
rather weak and in no sample yet, apparently, has one 
observed negative differential conductivity due to reso- 
nances associated with GaMnAs well layers. This hints 
at a significant concentration of defects, reminiscent of 
thin layers of amorphous Si where similar transport stud- 
ies have revealed size quantization effects but, to our 
knowledge, not negative differential conductivity. [23l [23| 
Experimental evidence indicating a coexistence of local- 
ized and extended Bloch-like states in bulk GaMnAs, in 
general, allows the prediction that, in thin layers, cer- 
tainly for < 3 nm, of GaMnAs extended states will be 
subjected to confinement effects (quantization and en- 
ergy shifts) while localized states will remain largely un- 
affected. This is similar to external magnetic-field effects 
on point defects or quantization effects in amorphous 
Si. [3, nil Assuming that no significant additional de- 
fects arise in GaMnAs heterostructures, this makes plau- 
sible experimental reports on quantum confinement ef- 
fects arising from (ferromagnetic) GaMnAs layers in thin 
heterostructures. [2TI. [22j Indeed, when one succeeds 
to incorporate high-quality magnetic layers in semi- 
conductor heterostructures strongly spin-dependent car- 
rier transmission can be predicted due to spin-selective 
tunneling. f26l] In magnetic resonant tunneling structures 
of high structural quality this spin splitting may be used 
for a realization of spin valves, spin filtering, and spin 
switching devices [21], [27l - [33j , all representing important 
ingredients for spintronic-based device technology. 

In several experiments ferromagnetism has been gener- 
ated in bulk GaMnAs by, electrically or optically, tailor- 
ing the hole density, (s^, [sl] In 2d-confined systems con- 
taining layers of Gai_a;Mna;As the magnetic order de- 
pends strongly on the local spin density, which can be 
influenced by the tunneling current, resulting in a bias- 
dependent exchange splitting. [H, 113 A spin-density de- 
pendent exchange splitting in ferromagnetic structures 
enriches the dynamic complexity by offering a mecha- 
nism for external electrical control of the ferromagnetic 
state. This is in contrast to structures comprising param- 
agnetic DMS, such as ZnMnSe, in which a giant Zeeman 
splitting of the bands is induced by applying an external 
magnetic field of the order of a few Tesla. 

Already nonmagnetic multi-well heterostructures ex- 
hibit interesting dynamic nonlinear effects which are 
based, however, on different physical mechanisms, such 
as the formation of electric field domains and the motion 
of charge dipoles through the structure. (38l - l4l| Recently 



it has been predicted that, in heterostructures containing 
paramagnetic DMS wells, this kind of phenomena can be 
controlled by an external magnetic field. [42|-|4^ Using an 
incoherent, sequential tunneling model we have proposed 
that ferromagnetic multi-well structures can generate ac 
spin currents, a phenomenon which originates from time- 
dependent inversion of the spin population in adjacent 
wells, [m 

In this article we investigate spin-selective hole trans- 
port in GaAs/AlGaAs/GaMnAs heterostructures within 
the limit of moderately thin samples with predominantly 
coherent transport characteristics. We apply a non- 
equilibrium Green's function formalism based on a tight- 
binding Hamiltonian for the electronic structure, includ- 
ing self-consistency regarding the charge density and 
the exchange splitting of the effective potential, as well 
as charge transfer to the contacts. Both the carriers' 
Coulomb interaction and the exchange coupling with the 
magnetic ions are described within a mean-field picture. 
Details of our model are exposed in Sect. |lll The mech- 
anism of electric control of magnetization switching is 
explored for two generic structures containing, respec- 
tively, one and two layers of Gai_2:Mn2:As. Results are 
given in Sect. IIIII We also provide a qualitative explana- 
tion for the occurrence of spin-polarized current oscilla- 
tions, predicted in an earlier paper (45j . and investigate 
the influence of spin flip processes at the interfaces on the 
total current spin polarization. Since disorder seems to 
play a major role in actual samples we study the effect of 
substitutional disorder on a qualitative level and discuss 
the robustness of the effects predicted here. Relevance 
to experiment is discussed. In particular, we can give an 
explanation for the absence of exchange splitting (mag- 
netization) under resonance bias condition reported in 
a recent experiment and identify characteristic features 
which may be explored in future experiments. Summary 
and conclusions are given in Sect. IIVI 

II. PHYSICAL MODEL 

The magnetic semiconductor heterostructure is de- 
scribed by a two-band tight-binding Hamiltonian for the 
heavy holes (J3 — ±3/2). It is given in the form 

+ '^t,^„^,\i,a){i + l,a'\+h.c., (1) 

where Ei^^ is the spin-dependent (a =^,\,= ±1) onsite 
energy at lattice site i, ti^cra' denotes the hopping-matrix 
between neighboring lattice sites, and h.c. abbreviates 
the Hermitian conjugate term. Spin conserving hop- 
ping gives a diagonal matrix ti^o-.o-' = t^aa', whereas 
spin flip processes can be taken into account by intro- 
ducing off-diagonal elements. The hopping parameter 
t = —lT?/{2m*a?') depends on the effective mass m* and 



3 



the lattice spacing a between to neighboring lattice sites. 
The onsite energy 



2 



(2) 



includes the intrinsic hole band profile Ui due to the band 
offset between different materials, the electrostatic po- 
tential (f) with e denoting the elementary charge, and the 
local exchange splitting A^. Near the band-edges this 
model is equivalent to an effective-mass model, however, 
it has the advantage that structural imperfections and 
spin-flip processes can readily be incorporated. More- 
over, it can be extended to arbitrary sophistication by 
introducing a larger set of basis functions. [26l. I46l - l48j 

Within a mean-field approach the exchange coupling 
between holes and magnetic impurities can be described 
by two interrelated effective magnetic fields, respectively, 
originating from a nonvanishing mean spin polarization 
of the ions' d-electrons (S^) and from the hole spin den- 
sity (sz) — (fit ~ '^4- )/2 -[IE 113,11^ The exchange splitting 
of the hole bands is then given by 



A(z) = -Jpdnimp{z){Sz){z) , 



(3) 



with z being the longitudinal (growth) direction of the 
structure, Jpd > is the coupling strength between the 
impurity spin and the carrier spin density (in case of 
GaMnAs p-like holes couple to the d-like impurity elec- 
trons), and rtimp(-2) is the impurity density profile of mag- 
netically active ions. Since the magnetic order between 
the impurities is mediated by the holes, the effective im- 
purity spin polarization depends on the mean hole spin 
polarization via 
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where ks denotes Boltzmann's constant, T is the lattice 
temperature, and Bs is the Brillouin function of order S, 
here with S ~ 5/2 for the Mn impurity spin. Combining 
Eq. ([3]) and Eq. (jH leads to a self-consistent effective 
Hamiltonian for the holes Hcs = —crA{z)/2 with 



A(z) = Jpdnimpiz)SBs 



SJpd[n-f{z) ~ nijz)] 
2kBT 



(5) 



Note that in thermodynamic equilibrium of a quasi 2D- 
sy stems, such as a quantum well, the hole spin density 
polarization (s^) is the key figure of merit for the appear- 
ance of ferromagnetism. 

Within a Hartree mean-field picture space-charge ef- 
fects are taken into account self-consistently by calculat- 
ing the electric potential from the Poisson equation. 



'^-e-^(l) = e[Naiz)-niz)], 



where e denotes the dielectric constant and Na is the 
acceptor density. The local hole density at site \i) is 
given by 



n(i) 
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with A being the in-plane cross sectional area of the 
structure, and fc|| denotes the in-plane momentum. The 
non-equilibrium "lesser" Green's function is calcu- 
lated from the equation of motion 

G< = G"^<G^ (8) 

where and G^ = [G^]"*" denotes the retarded and 
advanced Green's function, respectively. The scattering 
function S< — T,f + S< describes the inflow of particles 
from the left (l) and right (r) reservoirs [13] 



(9) 



where fo{x) = [1 + exp{x/kBT)]~^ is the Fermi dis- 
tribution function and ^/ and /i^, respectively, denote 
the quasi-Fermi energies in the contacts. The retarded 
and advanced self-energy terms S 



and 



= [S^]^ account for the coupling of the system re- 
gion to the left and right semi-infinite chains, for which 
an analytic expression can be derived. [sol [5l| The re- 
tarded Green's function is then given by 



G^ ^ [E + iTj- Hs- S^] ' 



(10) 



with irj being a positive infinitesimal imaginary part of 
the energy. 

Together with adjusting the Fermi energies relative to 
the band edges in the leads to ensure asymptotic charge 
neutrality '52i|, the band splitting given by Eq. ([5]), the 
Poisson equation Eq. (|6l7p . and the kinetic equations 
Eqs. ([5]) and (|10p have to be solved self-consistently until 
convergence to a steady-state solution is reached. Non- 
linearities in both Hartree and exchange term can give 
rise to multi-stable behavior, as will be discussed below. 
If this selfconsistency loop terminates with ferromagnetic 
ordering in the system, the effective one-particle poten- 
tial is different for spin-up and spin-down holes, thus 
leading to spin filtering in transmission. 

After obtaining the self-consistent potential profile the 
spin-dependent transmission probability T„'cr{E) from 
the left to the right reservoir, as a matrix element of 
the structure's S-matrix, can be calculated from special 
matrix elements of the retarded Green's function [481 



T , — T 



Jr,a'\G^{E-ra' M)? 
vi,,\G°{E-la, HI' 



(11) 



with Go denoting the free Green's function of the asymp- 
totic region, and vi^a and Vr.a, respectively, are the spin- 
dependent group velocities in the leads. G^{E;ra' ,la) 
is computed most conveniently by adding one layer after 
another which requires merely 2x2 matrix inversions for 
the present two-band model. [51| 

Finally, the steady-state current is obtained from scat- 
tering theory (generalized Tsu-Esaki formula). 



em*kBT 



AET,,ag{E) 



g{E) = In 



(27r)2ft3 

1 + exp [{in - E)/kBT] 



l-hexp[(/i^ -£;)/fcBT] 



(12) 
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The applied bias V = {m — /Zr)/e is determined by the 
difference in quasi-Fermi levels of the contacts. 

We would like to point out that we conduct a genuine 
non-equilibrium study whereby the quasi-Fermi level po- 
sitions are associated with the contacts. Self-consistency 
then leads to an effective, in general, spin-dependent 
one-particle potential. Thus one is not confronted with 
the question where to place the Fermi level in the GaM- 
nAs layers. Essential to confinement effects is the exis- 
tence of states near the top of the valence band edge of 
GaMnAs which have a coherence length of at least the 
layer thickness. Highly localized states, whether sepa- 
rated from or attached to the top valence band edge, will 
not be very sensitive to finite layer width. While in the 
bulk and thermal equilibrium the itinerant hole exchange 
model firmly relates hole density to Tc and the Fermi en- 
ergy, in a non-equilibrium tunneling situation this is dif- 
ferent. The key question is whether tunneling can induce 
a net hole spin polarization or not. As is shown below, 
we find that this depends on structural properties as well 
as on the applied bias. 
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III. 



RESULTS 



We start with a symmetric double-barrier structure 
containing a single GaMnAs quantum well and investi- 
gate the role of resonant hole tunneling on the magnetic 
state of the device. For the simulation we use generic pa- 
rameters for GaMnAs and GaAs: m* = 0.4mo, — 12.9, 
T4ar = 400 meV, /i, = /i^ = 80 meV, d ^ 20A, w = 25l, 
nimp = lxl020cm-3, = O.lSeVnm^ 53], r = 4.2K, 
where mg denotes the free electron mass, is the relative 
permittivity, Vbar is the bare barrier height of AlGaAs 
relative to GaAs, d and w, respectively, are the barrier 
and quantum well width. The thermal equilibrium po- 
sition of the Fermi energies /i; = fir was deliberately 
chosen close to the first resonance to promote ferromag- 
netic ordering in the well region at zero bias. The back- 
ground charge Na is assumed to be only about 10% of 
the Mn doping ni^p since GaMnAs is a heavily compen- 
sated system, most likely due to Mn interstitial or antisite 
defects. [^,|5jj The hole densities in the quantum well can 
be adjusted by the Mn doping level and the quasi-Fermi 
levels in the contacts. As can be seen from Eq. ([5]) the 
exchange splitting increases with the hole density in the 
case of a steady particle spin polarization. The value 
of the exchange coupling constant varies in literature to 
some extend Jpd ~ 0.04 — 0.16 eV nm. Since we use an 
optimistic value for Jpd, we assume only moderate MuQa 
doping in the well. Higher Mnca densities and smaller 
values of Jpd will give very similar results. 

Disorder effects in the GaMnAs layers are modeled by 
performing a configurational average over structures with 
randomly selected onsite and hopping matrix elements of 
the tight-binding Hamiltonian in the Mn doped region. 
For each specific Hamiltonian the transport problem (I- 
V curve) is solved self-consistently. The final result is 
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FIG. 1: (Color online) Spin-dependent transmission proba- 
bility of the double barrier structure at zero bias with and 
without disorder [at = 5%). 



obtained by averaging over all configurations. Typically 
300 configurations are used for one I-V curve. For the 
numerical simulation we assume a fixed 5% Mn concen- 
tration in the well and model substitutional disorder. If 
a Mn ion is present at a given lattice site in the well 
the onsite energy is shifted according to a Gaussian dis- 
tribution around a mean onsite energy-shift of 40 meV 
and a standard deviation of 20 meV, which are reason- 
able values according to experimental results which indi- 
cate either an impurity band slightly above the valence 
band edge or a defect-induced valence band tail. 
The nearest-neighbor hopping matrix elements for such 
a site are sampled according to a Gaussian between 5% 
and 25% standard deviation {(Tt) of its bulk value t. This 
model for substitutional disorder leads to a hybridization 
of quantum confined hole states, associated with bulk like 
valence band states, and localized defect states arising 
from Mn(3a sites. The degree of hybridization depends 
on layer thickness since it controls the position of the 
quantized heavy hole band relative to the energy of the 
localized Mn acceptor levels. This hybridization and the 
experimentally found Mn acceptor radius of about 2 nm 
calls for rather thin GaMnAs layers to ensure quantiza- 
tion effects in transport. [16] 

The calculated spin-filtering effect via distinct tunnel- 
ing probabilities for spin-up and spin-down holes arising 
from the exchange term is displayed in Fig. [T] in which 
the transmission probability at zero bias is plotted versus 
energy of incidence E. This figure also gives a qualita- 
tive account of the density of states in the GaMnAs well 
region discussed above. For an idealized GaMnAs layer 
which, at the valence band edge, is modeled as a GaAs 
layer plus exchange term, one obtains sharp spin dou- 
blets which are exchange-split by about 25-30 meV (see 
dashed versus solid lines in Fig. [T|). The state of zero 
spin polarization of holes represents an unstable equi- 
librium since, below Tc, the slightest perturbation in 
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spin-polarization drives the system into a partially or- 
dered lower energy state (spontaneous symmetry break- 
ing) due to the exchange interaction. The latter, in turn, 
accounts for different effective barrier profiles for spin- 
up and spin-down holes. It is this nonlinear effect that 
can be utilized to control the hole spin polarization and 
thus the favorable Mn spin orientation by structural de- 
sign and applied bias. Placing the Fermi level near the 
first heavy-hole resonance promotes this effect, similar 
to the formation of Cooper pairs near the Fermi edge of 
an interacting electron gas. 

Spin-selective tunneling into and out of the Mn doped 
wells, regardless of whether sequential or resonant, pro- 
motes hole spin polarization and, thus, alignment of the 
Mn spins as long as spin-depolarizing processes in the 
heterostructure are slow compared to the effective tun- 
neling rates. Furthermore, disorder which leads to spec- 
tral broadening of the resonances may suppress spin- 
selective tunneling. Inspection of Fig. [T] shows an asym- 
metric broadening and significant overlap of the trans- 
mission peaks under substitutional disorder, modeled as 
discussed above, which is particularly pronounced for the 
first heavy-hole resonance since it is most sensitive to po- 
tential fluctuations. The asymmetric ("anti-bonding") 
shift towards higher energies is due to hybridization with 
Mn acceptor levels below the conduction band edge. The 
latter do not contribute to resonant transport. Even at 
the moderate disorder for the effective hopping matrix 
element of 5 percent, a significant overlap in spin-up 
and spin-down resonance is obtained. Increased disor- 
der and/or spin-flip scattering will eventually wash out 
spin-selectivity in transmission and a destruction of fer- 
romagnetic ordering under bias must be expected since 
unpolarized holes are steadily fed into the GaMnAs re- 
gions. Exchange splitting at zero bias for 5% and 25%, 
respectively, is reduced to 33 meV and 23 meV. 

Clearly, our effective one-dimensional modeling of 
(substitutional) disorder must be viewed as a limited 
estimate since it corresponds to a cross-sectional av- 
erage of transport though uncorrelated effective linear 
chains. Correlations from disorder parallel to the het- 
erointerface will play a role in the establishing of coher- 
ence and ferromagnetic order in real structures relative 
to the idealized homogeneous mean-field model adopted 
here, since both ferromagnetic order and disorder effects 
are highly dependent upon spatial dimensionality. (55l . [56| 
Additional types of disorder from Mn clustering, Mn in- 
terstitials, etc. may be present in real structures. The 
role of disorder in the formation of ferromagnetic or- 
der in diluted magnetic semiconductors has been ex- 
plored theoretically and, remarkably, certain form of 
disorder has been predicted to promote ferromagnetic 
ordering. (Tsi [TqI In experiment, STM studies have given 
information on the nature of defects near the surface of 
GaMnAs samples. [lllist 

The current- voltage I-V characteristics, plotted in 
Fig. m reveals the typical hysteretic behavior of reso- 
nant tunneling diodes for an up- and down-sweep of the 
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FIG. 2: (Color online) /V'-characteristics of a magnetic double 
barrier structure due to heavy-hole associated bands with and 
without disorder. The solid (dashed) lines refer to a voltage 
up (down)-sweep without disorder (black lines) and moderate 
disorder of at — 5% (blue lines). The Jl^-curve is flattened 
at higher defect concentrations as indicated by the red dotted 
line at = 25% assuming a voltage up-sweep. As shown in the 
inset in the voltage range of = 0.07 — 0.09 V for the case of 
a high-quality sample no steady state is reached, suggesting 
the occurrence of dynamic effects. 




0.05 0.1 0.15 0.2 0.25 



V(V) 

FIG. 3: (Color online) The configuration averaged spin split- 
ting |A| in the quantum well as a function of the applied bias 
with and without disorder. 



applied bias. This well known intrinsic bistability oc- 
curs due to different charging of the well depending on 
the bias-sweep direction. Since our model ignores contri- 
butions from the light-hole band, associated resonances 
are missing in the plot. The latter are important due to 
in-plane non-parabolicity effects in narrow layers, how- 
ever, low-lying resonances associated with heavy and 
light holes generally are clearly separated in energy. For 
the present structure a light-hole-band resonance would 
be expected between the first two heavy-hole-associated 
resonances, thus strongly reducing the peak-to valley ra- 
tio and contributing spin ±1/2 holes to the Mn-doped 
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FIG. 4: (Color online) Averaged current spin polarization \Pj \ 
versus applied bias V with and without disorder. 



layers. It should be observed that only single resonance 
peaks are observed in the I-V characteristics in spite of 
spin doubles in the (zero-bias) transmission spectra. Fur- 
thermore, the drop in current beyond the first heavy-hole 
peak value (see insert in Fig. [2]), unlike in ballistic models 
for nonmagnetic tunneling structures (see second peak), 
is gradual even in the absence of disorder. This broaden- 
ing of the resonance can be attributed to ferromagnetic 
ordering avifay from the first resonance peak which tends 
to widen the bias window for meeting the resonance con- 
dition. 

Disorder effects diminish the peak-to valley ratio but 
regions of negative differential resistance are maintained 
for weak disorder. Since experiments have not shown 
negative differential conductivity in such a structure we 
have increasing disorder and find its disappearance at 
relatively high hopping disorder of about Ct = 25 % (see 
Fig. [2]). This indicates that defects other than Mn accep- 
tors are present in real samples. Further numerical stud- 
ies regarding this issue will be published elsewhere. [57| 

In Fig. [3] the average exchange band spin splitting |A| 
in the quantum well, characterizing its magnetic state, is 
plotted versus applied bias. It shows that ferromagnetism 
can be controlled by the applied bias in this structure 
near the first current peak, remarkably, even when disor- 
der is sufficiently strong to suppress negative differential 
conductivity. At zero bias ferromagnetic ordering is ener- 
getically preferred since the Fermi [ii = [ir level is located 
close to the edge of the first heavy-hole subband. As the 
bias is increased tunneling into the upper doublet state 
becomes allowed from the emitter side reducing the net 
hole-spin polarization (in spite of increasing hole den- 
sity) and the effective exchange field decreases to zero 
and both spin-up and spin-down subbands go into reso- 
nance. Note that under moderate bias both emitter and 
collector contribute to the population of the well region. 
As the bias increases resonant population from the emit- 
ter gets shut off and hole polarization is determined by 
the collector leading once more to a build-up of the ex- 



change field for a bias regime between 0.08 V and 0.2 V, 
when finally the collector quasi-Fermi level drops below 
the hole subbands and the well region becomes almost 
depleted of holes. For higher bias no further sponta- 
neous magnetization has been obtained within our self- 
consistency loop. The overall feature of the bias depen- 
dence of the exchange splitting thus somewhat resembles 
its behavior versus temperature, with "T = Tc" corre- 
sponding to a bias of about 0.18 V. It arises from the 
fact that it is the number of spin-polarized holes which 
determines the maximum spontaneous magnetization for 
given Mnca concentration. A simple model for the de- 
pendence of the Curie temperature in resonant tunnel- 
ing systems has been given by one of us before, [ssj The 
voltage-dependence of the Curie temperature under res- 
onant tunneling has also been studied before. [H^ 

The displayed build-up and destruction of ferromag- 
netic order as a function of applied bias can be further 
understood by the exchange interaction which is medi- 
ated by spin-polarized holes. In an ideal 2D particle sys- 
tem with parabolic dispersion there is no energy gain by 
magnetic ordering due to the constant density of states 
associated with each spin subband: energy gained by 
lowering one subband is exactly cancelled by raising the 
other. However, here we deal with a 3D heterostructure 
which favors a spin ordered state when the quasi-Fermi 
level lies near (within about half of the maximal exchange 
splitting) the bottom of a well subband resonance. If the 
temperature in the contacts is sufficiently low, one sub- 
band after the other will go through resonance. Thus, 
when only the lower spin-subband is in resonance holes 
in the magnetic well will tend to be be spin-polarized. 
However, as bias is increased eventually the subband with 
opposite spin orientation will also go into resonance thus 
reducing spin polarization and magnetic ordering in the 
GaMnAs layer. When, for a given bias, the well region 
cannot be populated (lack of hole density of states) or no 
energy gain can be drawn from ferromagnetic ordering, 
loss of the latter will result. 

Interestingly, in the voltage range oiV = 0.07 — 0.09 V 
no steady state solution can be found for the low disorder 
sample case. Instead the solution for the magnetization 
is oscillating, as shown in Fig. [3l suggesting the occur- 
rence of dynamic effects. This behavior can be under- 
stood qualitatively as follows: Figure [5] shows a contour 
plot of the local density of states, for an applied bias 
V — 0.085 V lying in the critical voltage range. The 
self-consistent band profile is indicated by the solid line. 
For the emitter Fermi energy of /i/ = 0.08 eV only the 
two ground state (potentially spin-split) subbands in the 
quantum well participate in the tunneling transport. At 
this bias condition and hole spin polarization the low- 
est (spin-up) subband may be populated by holes from 
the collector side, whereas the spin-down level is almost 
empty since it cannot be reached elastically by either 
emitter or collector. Since the (steady-state) band split- 
ting A is proportional to the spin polarization (n-j- — n\J) 
the well magnetization increases with spin polarization. 
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FIG. 5: (Color online) Logarithmic local density of states 
(LDOS) as a function of energy at the bias V — 0.085 V. The 
self-consistent band profile is indicated by the solid line. The 
spin-splitting of the quasi-bound states is clearly visible. 
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FIG. 6: (Color online) Current spin polarization Pj versus 
applied bias V taking into account spin flips at the hetero- 
interfaces. The polarization is diminished for an increasing 
spin flip probability p becoming unpolarized for p = 1/2. 



pushing the spin-down level upwards in energy. At some 
point holes can start to tunnel from the emitter side into 
the spin-down level. This in turn decreases the total 
spin polarization and, hence, effectively pushes the spin- 
down level back below the emitter's band edge. From 
there, the process starts anew, leading to an oscillatory 
behavior in well magnetization, tunneling current, and 
spin polarization, [isj 

Although the I-V curves in Fig. [5] do not display spin- 
split resonance peaks, but merely a broadening of the 
resonance, the steady-state current at low bias is spin 
polarized as shown in Fig. 2] As the bias is increased 
from zero, current spin polarization is reduced and re- 
versed before it drops to zero through resonance. Above 
resonance current spin polarization reemerges (due to the 
action of the collector) and once more changes sign be- 
fore dropping and remaining at zero in one-to-one agree- 



ment with the behavior of the exchange field. Although 
resonance peaks in the I-V curve my be suppressed by 
disorder, see Fig. [21 the bias-dependence of spin polar- 
ization in the current may persist and may be observed 
in experiment as a bias-dependent spin valve. 

In order to study qualitatively the influence of spin flip 
processes at the hetero-interfaces on the total current 
spin polarization at the collector side, Pj = (j^-|- + j^i — 
jit-ju)/j with j = J2a<7' we introduce off-diagonal 
hopping matrices Vi^aa' in the tight-binding Hamiltonian. 
In general for N interfaces there are 2^ different flip con- 
figurations. For each of them a simulation is performed 
and the results are finally averaged by weighting with the 
probability for the occurrence of the configuration. In the 
case of a double-barrier structure we have four hetero- 
interfaces, giving 16 configurations. However, flipping at 
the first barrier interfaces is inefficient, since it does not 
change the total current or spin polarization. Single flip- 
ping at the third or fourth interface does also not modify 
the total current density but inverts the spin polariza- 
tion to —Pj- By introducing single spin flip probabilities 
Pi, {i — f , . . . , N) at the interface i, the probability of a 
flipping process at the second barrier is then given by 
Pflip = P3(l — Pi) + (1 — P3)P4- Hence, the mean spin 
polarization results in 



(13) 



The bias-dependent current spin polarization for different 
spin flip probabilities (assuming ps — P4: — p) is plotted 
in Fig. ini The spin polarization decreases for increasing p 
with {Pj)[p] = {Pj)[l~p] reaching its minimum {Pj) = 
for p = 1/2. From this analysis we conclude that our 
results will not be altered significantly when the spin- 
orbit interaction is taken into account in the analysis. 

While spin-selective hole tunneling may allow electric 
control of ferromagnetic order, tunneling spectroscopy, 
in turn, provides a sensitive experimental tool for ex- 
ploring the electronic structure of mesoscopic semi- 
conductor systems. [sO] Recently, tunneling spectroscopy 
experiments have been performed on thin layers of 
GaMnAs. [l^ The authors have verified ferromagnetic or- 
dering in their samples (with Mn concentration of typi- 
cally X w5 to 15 % and layer thickness ranging from 4 to 
20 nm) and have measured their respective Curie temper- 
ature. Their measurements indicate that Mn induced de- 
fect states remain separated from the GaAs-like valence 
band edge as evidenced by a pinning of the Fermi level. 
Furthermore, they find clear signatures of quantization 
effects in the transmission spectra of their samples and re- 
port an absence of spin-splitting in the resonances which 
they can fit to a GaAs-like k.p model, including light-hole 
states. We believe that these experimental findings com- 
pare favorably with the general features of our results. 
Moreover, we have provided an explanation for the ob- 
served absence of ferromagnetic ordering near resonance 
in spite of ferromagnetic behavior of the sample at zero 
bias. It would be interesting to perform spin-sensitive 
tunneling spectroscopy on these samples since, accord- 
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FIG. 7: IV-characteristics of a three-barrier structure with 
two coupled quantum wells made of GaMnAs. At the 
current maxima resonance conditions are fulfilled, i.e., the 
quasi-bound states of the adjacent wells become energetically 
aligned. The inset shows the local density of states at the ap- 
plied bias V — 0.03 corresponding to the first current maxima. 
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FIG. 8: (Color online) Maximum exchange splitting A max m 
the first (solid) and second well (dashed line) as a function of 
the applied bias. 



ing to Fig. ini such a measurement gives more detailed 
information about the bias dependence of ferromagnetic 
ordering than the I-V curve and its derivatives. This can 
test the prediction that ferromagnetic order which can be 
achieved at zero bias is destroyed near resonance would 
be verified, and that electric switching back and forth 
between the ferromagnetic and paramagnetic state can 
be achieved. 

We now explore feasibility of selective magnetization 
switching among several magnetic layers of high struc- 
tural quality. We investigate a three-barrier structure 
with two adjacent GaMnAs quantum wells, choosing an 
asymmetric structure with the second well being thinner 
than the first one {wi — 25 A, W2 = 20 A). All other 
parameters are as in the previous structure. Quantum 



confinement gives rise to a higher ground state energy 
in the second quantum well at zero bias. The resonant 
alignment of the ground state subbands of the two wells 
is therefore achieved at a finite voltage as shown in the 
inset of Fig. [3 corresponding to the first maximum in 
the current-voltage characteristics at about V = 0.02 V 
which is plotted in Fig. [T] The second current maximum 
result from the resonance of the first excited state sub- 
bands of both wells. Next to possible exchange splitting 
the finite separating barriers cause the energy levels in 
the two quantum wells to further split into bonding and 
antibonding subband states. However, for our structure 
the middle barrier is too thick and the natural energy 
broadening of the quasi-bound states is too large for re- 
solving this additional splitting in the local density of 
states. 

Having two coupled quantum wells allows one to realize 
several magnetic configurations. The maximum (steady- 
state) exchange splitting of the two wells as a function 
of the applied bias is plotted in Fig. [51 revealing three 
different regions. For low voltages both wells are mag- 
netized due to the build-up of spin polarization in the 
wells due to resonance of the ground state subband lev- 
els with populated reservoir states. Exchange causes a 
relative shift in the density of states for spin-up and spin- 
down holes which, in turn, stabilizes ferromagnetism in 
both layers. When the second well goes off resonance at 
around V > 0.03 V the accumulated spin polarization 
in the second well is preserved, since for voltages up to 
y « 0.1 V the collector Fermi energy is still higher 
than the bottom of its ground state subbands thus main- 
taining spin polarization. For voltages in the interval 
0.12 V< V < 0.33 V the first well remains magnetic, 
whereas the second well becomes nonmagnetic, since the 
ground state subbands are no longer filled from the col- 
lector side. At sufficiently high bias, V > 0.33 V, also the 
first well becomes demagnetized since holes can no longer 
resonantly populate its two lowest (now degenerate) sub- 
bands from either emitter or collector thus resulting in a 
completely nonmagnetic structure. 

Several simplifying assumptions have been made in 
the present analysis which should, just as well as ex- 
perimental aspects, be discussed. The present model is 
based on an effective-mass-like two-band approach for 
the heavy holes in the structure. This approximation 
should at least qualitatively be correct since the applied 
bias is kept below typically 0.2 V and most of the phe- 
nomena discussed here occur at lower bias. Thus it can 
be expected, that this model describes effects qualita- 
tively correct. We are currently working on more real- 
istic tight-binding formulations using a significantly in- 
creased number of basis states in conjunction with den- 
sity functional plus dynamic mean-field models to arrive 
at a more detailed and realistic electronic structure. [Slf 
Impurity scattering effects have been accounted for on a 
phenomenological level within the TB model. Our bal- 
listic model neglects electron-phonon scattering within 
the heterostructure altogether and the electron-electron 
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interaction is described within mean-field theory. In thin 
structures, such as the ones studied here where effective 
tunneling rates are higher than carrier-phonon scattering 
rates and optical phonon transitions are suppressed en- 
ergetically the former assumption should be rather well 
fulfilled and not alter significantly subband population 
within the heterostructure. Electron-electron scattering 
may play role, however, as long as it does not involve 
spin-flip processes should not influence our basic conclu- 
sions much. 

Clearly the effects studied here require low tempera- 
tures, for one to favor ferromagnetic ordering and, sec- 
ondly, to preserve strong hole-spin polarization in the 
carrier injection process. It is well known that, at least 
at low temperatures, structural imperfections are the 
main source for reduction of nonlinear effects^such as the 
peak-to-valley-ratio in the IV curve . p7l |63 - |64| It is most 
likely the difhculty in clean sample preparation which 
has slowed experimental progress on thin-layer semimag- 
netic semiconductor heterostructures. High quality dop- 
ing profiles and high quality interfaces must be achieved 
within one growth process. [U [3l| Growth of good 
quality DMS layers needs low temperature molecular 
beam epitaxy which, however, adversely affects interface 
quality. Usually thin GaAs spacer layers are inserted 
to smooth the surfaces. [2l|, jsj Furthermore, GaMnAs 
layers must be thick enough to support ferromagnetism. 
Qualitatively, all structural imperfections lead to broad- 
ening of resonances. Once the latter becomes comparable 
to the (theoretical) maximum of the exchange energy in- 
duced spin-splitting, spin-selective tunneling and, hence, 
tunneling-induced control of magnetic ordering may be 
suppressed. Even in the presence of disorder, as long as 
it does not go hand in hand with strong spin-hip pro- 
cesses, achieving bias-control of hole-spin polarization 
in the GaMnAs layers should allow one to manipulate 
magnetization. 

IV. CONCLUSIONS AND OUTLOOK 

In summary, we have used a ballistic steady-state 
transport model to investigate bias-induced magnetic 
multi-stability in AlGaAs/MnGaAs quantum well struc- 
tures. Ferromagnetic exchange, as well as the hole 
Coulomb interaction are treated within self-consistent 
mean-field approximation. Substitutional disorder 
is treated phenomenologically within a tight-binding 
model. 

Our studies indicate that in these systems ferromag- 
netic ordering can be controlled selectively by an exter- 
nally applied bias. The underlying mechanism is found 
in spin-selective tunneling due to the anti-ferromagnetic 
exchange interaction between itinerant heavy holes and 
localized Mn d-electrons. In structurally suitably de- 
signed heterostructures the applied electric bias allows 
control of the ferromagnetic state, as well as electric and 
spin current density. 



In the simplest structural case in form of a double bar- 
rier structure containing a GaMnAs well, we predict that 
ferromagnetic ordering in the well, when present at zero 
bias, is lost under bias near the first heavy-hole reso- 
nance, allowing a switching back and forth between the 
magnetic and nonmagnetic state in the well. In GaM- 
nAs multi-well structures we predict that the loss of fer- 
romagnetic order can be engineered structurally to occur 
at different applied bias for the individual layers. 

Within our model we are able to provide a possible 
explanation for the absence of exchange splitting near 
resonances, as observed in recent tunneling spectroscopy 
measurements on thin GaMnAs layers. [1^1 We generally 
predict that ferromagnetic order which may be achieved 
in GaMnAs quantum well layers under zero bias tends to 
be destroyed under resonance condition since the well re- 
gion then is swept by unpolarized holes. Under favorable 
conditions detailed in the main text, ferromagnetic order 
may be reestablished above resonance. Such a behav- 
ior should be revealed experimentally by spin-sensitive 
tunneling spectroscopy. [65| 

In previous work based on a complementary time- 
dependent sequential tunneling model including intra- 
well scattering we have predicted that, under specific 
bias conditions, the interplay of transport and magnetic 
properties can result in robust self-sustained charge and 
magnetization oscillations. [45| The present model, albeit 
based on the resonant-tunneling picture, backs the pos- 
sibility of such phenomena by predicting bias regions in 
which no steady-state solution for the current exists. 

Disorder and spin-flip effects have been modeled on 
a phenomenological level. We find that disorder due to 
Mn taking a Ga site alone should not suffice to destroy 
spin-selective tunneling, nor should spin fiips at a rate 
expected in these structures, for example from the spin- 
orbit interaction, significantly suppress spin polarization 
of the steady-state current. As expected, our analysis 
does show that disorder and spin fiip processes do reduce 
the total average current spin polarization, however, not 
as efficiently as the resonance peaks in the I-V curve. 

We conclude that multi-well structures containing 
GaMnAs layers may allow one to realize various bias- 
dependent magnetic configurations. While the current 
investigation of bias induced effects considers only bias 
in longitudinal direction, i.e., a 2-terminal configuration, 
applying additional gates in transverse direction (multi- 
terminal configuration) should allow for an additional 
control knob to move spin-split subbands in and out of 
resonance with the contact states and/or to inject spin- 
polarized holes into the Mn-doped regions. Such a struc- 
ture has been studied in a recent experiment . [sH 
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